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ABSTRACT 


Quantitative  information  about  human  frequency  response 
functions  in  the  control  of  a  piloted  aircraft  would  provide  a 
rational  basis  for  stability  and  control  criteria.  Were  the 
pilot  a  linear,  time  invariant  system,  the  ratio  of  the  spectral 
densities  of  his  output  signal  to  his  input  signal  would  be  the 
square  of  the  transfer  function  describing  the  pilot.  In  addi¬ 
tion,  the  computation  of  these  spectral  densities  affords 
valuable  insights  into  the  response  process  under  study.  In 
an  effort  to  examine  the  invariance  of  human  system  parameters 
an  attempt  was  made  to  change  an  operator's  tracking  responses  by 
using  two  different  sets  of  instructions  for  the  tracking  of  the 
same  statistical  input.  Two  subjects  were  used  in  a  simple 
compensatory  position  tracking  device.  Although  the  amplitude 
response  characteristics  do  not  conclusively  prove  or  disprove 
the  invariance  of  the  response  patterns,  the  spectral  densities 
of  the  signals  provide  insight  into  the  manner  in  which  the 
type  of  instructions  affected  the  response.  The  value  of 
future  research  is  indicated. 
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SECTION  I 


INTRODUCTION 


The  purpose  of  this  report  is  to  illustrate  an  application  of 
the  spectral  density  approach  to  the  study  of  human  responses  in  a 
perceptual  motor  problem,  fi]  The  particular  aspect  of  this  problem 
with  which  this  report  will  deal  is  the  influence  of  instructions  on 
the  performance  of  a  tracking  task.  The  type  of  data  presentation 
resulting  from  the  computation  of  the  spectral  density  of  the  time 
functions  in  a  tracking  task  would  appear  to  be  appropriate  to  the 
determination  of  differences  in  motor  response  patterns  due  to  what 
are  presumed  to  be  differences  in  set  to  respond  associated  with 
different  instructions  to  the  subject.  Knowledge  of  the  spectrum 
of  an  input  to  a  device  under  study  and  the  resulting  output's 
spectrum  enable  conclusions  to  be  reached  about  the  type  of  opera¬ 
tion  which  the  device  performs  on  input  signals. 

The  raw  data  for  the  spectral  density  analysis  were  obtained 
by  means  of  a  compensatory  position  tracking  device,  the  control 
ratio  of  which  was  fixed  at  one  centimeter  of  deflection  on  the 
scope  face  per  degree  of  control  stick  motion.  This  device  was 
was  designed  to  test  the  applicability  of  a  spectral  density 
analysis  to  the  reduction  of  data  from  a  dynamic  aircraft  simu¬ 
lator.  Future  and  more  ambitious  experiments  involving  the 
methods  outlined  in  this  report  will  be  performed  with  the  dynamic 
simulator.  It  is  hoped  that  these  future  tests  will  serve  to  shed 
more  light  on  the  response  characteristics  of  pilots. 
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SECTION  II 


EXPERIMENTAL  DETAILS 


A.  APPARATUS  AND  INPUT  PROGRAM 

The  apparatus  used  in  this  experiment  was  a  one-dimensional  track¬ 
ing  device  as  illustrated  in  Figure  1.  This  device  was  designed  to 
serve  as  a  simplified  approximation  to  a  dynamic  aircraft  simulator 
for  the  purpose  of  providing  a  convenient  means  for  planning  research 
using  the  dynamic  simulator.  The  display  consisted  of  a  cathode-ray 
•scope  mounted  approximately  2S  inches  from  the  subject’s  eyes,  A 
pip  on  this  ’scope  executed  discrete  left  and  right  steps  in  a  fixed 
horizontal  line  on  the  'scope  face.  The  subject’s  task  was  to  center 
this  pip  on  a  vertical  fiducial  line.  The  joystick  control  permitted 
position  tracking  of  this  spot.  The  apparatus  controls  were  designed 
to  allow  the  programming  of  20  square  pulses  of  durations  which 
could  vary  in  eight  discrete  steps  from  0,25  to  1.91  seoonds.  A 
program  was  devised  so  that  the  pulse  durations  were  selected  from 
a  Poisson  distribution  of  zero  crossings,  whereas  the  pulse  ampli¬ 
tudes  were  alternately  to  the  left  and  right  with  mean  magnitudes 
of  one  centimeter  in  each  direction.  The  average  pulse  duration  was 
•75  seconds  and  the  total  length  of  the  20  pulse  program  was  a  little 
over  15  seconds.  The  individual  amplitudes  were  selected  from  a 
Gaussian  distribution  of  mean  unity  and  standard  deviation  0,25  om. 

Recordings  were  made  on  a  Brush  recorder  of  the  control  stick 
position,  ,  the  instantaneous  pip  position,  £  ,  and  the  step 
function  program,  Figure  2  illustrates  a  segment  of  a  typical 

record.  The  rounding-off  of  the  square  input  pulses  is  due  to  the 
filtering  action  of  the  recorder’s  simplifiers. 


5TICK,  ff0 


Figure  2,  Ten  Seconds  of  Recordings  of  £  ,  0;  ,  and  $.  for  Range  1. 
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SECTION  II 


EXPERIMENTAL  DETAILS 


A.  APPARATUS  AND  INPUT  PROGRAM 

The  apparatus  used  in  this  experiment  ms  a  ane-dinensional  track¬ 
ing  device  as  illustrated  in  Figure  1.  This  device  was  designed  to 
serve  as  a  simplified  approximation  to  a  dynamic  aircraft  simulator 
for  the  purpose  of  providing  a  convenient  means  for  planning  research 
using  the  dynamic  simulator.  The  display  consisted  of  a  cathode-ray 
'scope  mounted  approximately  22  inches  from  the  subject's  eyes#  A 
pip  on  this  'scope  executed  discrete  left  and  right  steps  in  a  fixed 
horizontal  line  on  the  'scope  face.  The  subject’s  task  was  to  center 
this  pip  on  a  vertical  fiducial  line.  The  joystick  control  permitted 
position  tracking  of  this  spot#  The  apparatus  controls  were  designed 
to  allow  the  programming  of  20  square  pulses  of  durations  which 
could  vary  in  eight  discrete  steps  from  0.25  to  1.91  seoonds#  A 
.program  was  devised  so  that  the  pulse  durations  were  selected  from 
a  Poisson  distribution  of  zero  crossings,  whereas  the  pulse  ampli¬ 
tudes  were  alternately  to  the  left  and  right  with  mean  magnitudes 
of  one  centimeter  in  eaoh  direction.  The  average  pulse  duration  was 
#75  seconds  and  the  total  length  of  the  20  pulse  program  was  a  little 
over  15  seconds.  The  individual  amplitudes  were  selected  from  a 
Gaussian  distribution  of  mean  unity  and  standard  deviation  0.25  cm. 

Recordings  were  made  on  a  Brush  recorder  of  the  control  stick 
position,  Sm  ,  the  instantaneous  pip  position,  £  ,  and  the  step 
function  program,  .  Figure  2  illustrates  a  segment  of  a  typical 
record.  The  rounding-off  of  the  square  input  pulses  is  due  to  the 
filtering  action  of  the  recorder’s  amplifiers. 


Figure  2.  Ten  Seconds  of  Recordings  of  £  ,  9;  ,  and  for  Range  1. 
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The  choice  of  the  foregoing  program  was  determined  by  a  desire 
to  present  a  visual  display  which  would  have  a  fairly  smooth  spectral 
density  and  which  would  provide  a  task  requiring  the  operator  to 
pay  close  attention  to  the  pip  on  the  scope  face.  Figure  3  is  a 
plot  of  the  spectral  density  of  the  selected  program  of  square  pulses. 


Figure  3.  Spectral  Density  for  the  Program,  Q j  . 


B.  EXPERIMENTAL  PROCEDURE 

Two  male  subjects  were  used  in  this  experiment,  subject  #1  was 
a  retired  military  man  who  had  at  least  thirty  years  of  flying 
experience.  Subject  #2  was  a  man  in  his  early  twenties  who  had  no 
flight  experience.  The  same  program  of  square  pulses  was  presented 
to  each  subject,  and  each  subject  was  given  intensive  training  in 
the  tracking  problem.  The  experimental  variable  was  the  subject's 
instructions.  On  the  first  day  of  the  experiment  subject  #1  was 
instructed  as  follows:  "It  is  of  crucial  importance  to  this  study 
that  you  pay  careful  attention  to  these  instructions.  Return  the 
pip  to  the  vertical  line  as  accurately  as  you  can."  On  the  same 
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day  subject  #2  was  given  instructions  which  differed  from  the  fore¬ 
going  only  in  that  the  word  "rapidly"  was  substituted  for  the  word 
"accurately".  Prior  to  recording  the  three  time  functions,  ®o>  ®i» 
£  for  each  subject,  the  subjects  received  a  familiarization  train¬ 
ing  of  a  few  minutes  with  the  simplified  simulator.  Following 
this  session,  the  subject  underwent  four  training  sessions  of  five 
minutes  each,  separated  by  intervals  of  about  one  hour.  During 
these  sessions  each  subject  repeatedly  tracked  the  same  program 
under  his  particular  instructions,  which  were  repeated  at  the 
onset  of  each  training  period.  An  hour  after  the  last  training 
session  the  subject  was  once  again  given  the  appropriate  instruc¬ 
tions  and  presented  a  succession  of  five  20-step  programs.  The 
three  measured  signals  were  recorded  for  the  last  two  of  these  pro¬ 
grams,  and  these  recorded  signals  were  considered  to  be  representa¬ 
tive  of  a  trained  subject. 

On  the  second  day  of  the  experiment  the  same  procedure  was 
adhered  to  except  that  subject  #1  was  instructed  to  track  for 
speed  and  subject  #2  was  instructed  to  track  for  accuracy. 


DATA  REDUCTION 

In  order  to  convert  the  time  function  records  to  spectral  den¬ 
sities,  it  was  necessary  to  obtain  their  autocorrelation  functions. 
This  process  was  accomplished  by  means  of  a  mechanical  analogue 
device  at  the  MIT  Servomechanisms  Laboratory.  [2^1  The  cosine 
transform  was  obtained  by  means  of  a  2h  harmonic  analyzer  at  the 
MIT  Instrumentation  Laboratory.  The  details  of  the  data  reduction 
procedure,  as  well  as  an  outline  of  the  theoretical  background  of 
the  procedure  is  to  be  found  in  AFTR-6723.  \_l3 
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SECTION  III 


RESULTS 


Figures  U  and  5  show  the  spectral  density  plots  for  0o, 
and  £  for  subject  #1  under  instructions  for  accuracy  and  instruc¬ 
tions  for  speed.  Figures  6  and  7  are  analogous  plots  for  subject 
#2.  The  confidence  limits  denoted  by  the  vertical  lines  represent 
pessimistic  approximations  to  spreads  of  20  on  either  side  of  the 
mean  due  to  errors  in  the  mechanical  determination  of  the  autocorre¬ 
lation  function,  and  do  not  include  the  errors  due  to  the  sampling 
and  truncation  inherent  in  the  spectral  density  computations.  The 
spectral  densities  have  been  smoothed  according  to  Tukey's  procedure. [1,30 

It  will  be  noticed  from  Figures  U  and  5  that  there  is  little 
difference  between  the  response  of  subject  #1  under  instructions 
for  accuracy  and  under  instructions  for  speed.  The  explanation  of 
this  fact  may  conceivably  be  based  on  the  following*  subject  #1  may 
have  interpreted  the  words  accuracy  and  speed  in  the  framework  of  his 
past  experience  with  tracking  devices,  and  he  may  have  realized  that 
in  order  to  achieve  a  minimum  error  he  would  have  to  track  as  rapidly 
as  possible. 

Subject  #2,  on  the  other  hand,  gives  evidence  of  a  different 
type  of  response  for  instructions  to  track  for  accuracy  as  compared 
with  instructions  to  track  for  speed.  This  can  be  seen  by  comparing 
Figure  6,  which  represents  tracking  for  accuracy,  with  Figure  7, 
which  represents  tracking  for  speed.  In  Figure  6  the  error  signal’s 
spectral  density  is  predominantly  in  the  low  frequency  range  with  a 
peak  at  about  3/1*  cycles  per  second}  thus  it  might  be  implied  by 
comparing  the  error  spectral  density  with  the  program's  spectral  den¬ 
sity  at  higher  frequencies  that  the  tracker  was  not  introducing  any 
of  the  higher  frequencies  into  the  error  signal.  Subject  #2  was,  in 
effect,  smoothing  his  response  and  lumping  his  output  in  the  lower 
frequencies  when  instructed  to  track  for  accuracy.  This  gives  insight 
into  subject  #2's  interpretation  of  the  word  "accuracy" .  On  the  other 
hand.  Figure  7,  which  is  subject  #2's  response  pattern  under  instruc¬ 
tions  for  speed,  demonstrates  that  when  instructed  to  respond  for  speed, 
subject  #2  effectively  eliminated  the  smoothing  of  his  response  and  re¬ 
distributed  his  output  over  higher  frequencies  with  somewhat  of  a  peak 
at  a  frequency  of  one  cycle  per  second.  Subject  #2  might  be  thought 
of  as  attempting  to  square  up  his  response  to  the  step  inputs  by 
introducing  more  high  frequency  components  in  his  output. 
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Figure  5-  Spectral  Densities  of  €  ,  ,  for  Subject  1,  Instructed  for  speed. 
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Figure  7.  Spectral  Densities  of  £  ,09,S)  for  subject  2,  Instructed  for  Speed. 


It  is  of  interest  to  compare  the  root-mean-square  error  for 
the  two  subjects  under  the  two  instructions.  Consider  the  follow¬ 
ing  table  of  root-mean-square  errors  in  centimeters  on  the  scope 
face. 


Subject  #1  Subject  #2 


Instructions  for  Accuracy  .53  .61* 


Instructions  for  Speed  .$U  .5U 


It  is  seen  that  subject  #2  actually  was  doing  considerably 
poorer  with  regard  to  his  root-aean-square  error  when  he  tracked 
for  accuracy  than  when  he  tracked  for  speed.  Perhaps  the  reason 
for  the  difference  between  his  performance  and  the  performance  of 
subject  #1  lies  in  the  order  in  which  the  experiments  were  performed. 
Subject  #2  was  instructed  for  speed  on  the  first  day  and  accuracy  on 
the  second.  Hence  it  is  conceivable  that  being  a  cooperative  and 
intelligent  subject,  he  realized  that  a  difference  was  suspected  to 
exist  between  his  performance  on  the  two  days  and  thus  he  attempted 
to  somehow  distinguish  his  responses  to  the  accuracy  instructions  on 
the  second  day  from  his  responses  to  the  speed  instructions  on  the 
first  day.  Subject  #2  was  in  a  sense  forced  to  interpret  accuracy 
as  meaning  something  distinguished  from  speed.  Subject  #1,  however, 
was  instructed  for  accuracy  first  and,  responding  as  he  did,  he 
could  not  produce  a  different  type  of  response  under  instructions 
for  speed. 

The  results  can  be  presented  in  different  form.  Since  the 
intent  of  this  study  is  to  determine  means  for  obtaining  a  transfer 
function,  it  would  be  useful  to  compute 

h(p)l  «  I  4..(«)  ,  (i) 


in  order  to  determine  what  regularities  this  function  preserves 
despite  changes  in  instructions,  as  well  as  to  see  whether  or  not 
this  function  has  a  shape  which  would  look  reasonable.  There  are 
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two  independent  ratios  which  can  be  oomputed,  but  the  computation  of 
three  ratios  in  all  serves  to  present  the  results  more  clearly.  Con¬ 
sider  Figure  0  which  is  a  simplification  of  the  traoking  problemj 


Figure  S.  Simplified  Diagram  of  Tracking  Problem. 


Oo/ 0}  is  the  amplitude  ratio  across  the  open  tracking  loop,  So/e 
is  the  amplitude  ratio  across  the  human  operator’s  section  of  the 
loop  and  is  the  amplitude  ratio  for  the  closed  loop. 


Figures  9»  10,  and  11  present  the  foregoing  ratios  for  subjects 
#1  and  #2  under  instructions  for  speed;  whereas  Figures  12,  1J,  and  11+ 
are  these  ratios  for  subjects  jfl  and  #2  under  instructions  for  accuracy. 
The  fiducial  limits  represent  2tf*  confidence  limits  based  on  the  con¬ 
fidence  limits  in  Figures  i+,5,6,  and  7»  [l]  The  bandwidths  over 

which  the  amplitude  ratios  were  plotted  were  limited  so  that  the  2.0" 
confidence  limits  did  not  beoome  so  large  as  to  render  the  plots 
meaningless. 


The  amplitude  ratios  present  certain  features  which  demand  ex¬ 
planation.  The  amplitude  ratio  |  Oo/&;\  in  Figures  9  and  12  would 
reasonably  be  expected  to  equal  unity  for  low  frequencies,  since 
one  would  expect  that  the  tracker  could  position  a  pip  perfectly 
under  d-c  conditions.  Figures  9  and  12  depart  from  unity  to  an  ex¬ 
tent  which  cannot  be  explained  in  terms  of  the  computational  error 
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confidence  bands.  If,  however,  one  considers  Tukey's  £3]  concept 
of  assigning  degrees  of  freedom  to  spectral  densities  and  thus 
arriving  at  confidence  limits  for  the  sampled  estimate  of  the  time 
spectrum,  it  is  possible  to  set  additional  confidence  bands  about 
the  amplitude  ratios.  The  spectral  densities  on  Figures  U,5,6  and 
7  represent  28  degrees  of  freedom  til  sund  for  this  many  degrees  of 
freedom  Tukey  has  shown  that  the  ratio  of  the  ratios  of  observed 
powers  in  two  frequency  bands  to  the  actual  ratio  of  powers  in 
these  two  bands  will  exceed  the  limit  of  1.88  ten  percent  of  the 
time.  Since  the  amplitude  ratios  are  the  square  root  of  the 
power  spectra  ratios,  the  90%  level  limit  of  the  ratio  of  the  ratios 
would  be  1.37. 

This  fact  makes  the  deviation  of  the  amplitude  ratio  from 
unity  at  low  frequencies  somewhat  less  unreasonable.  A  further 
explanation  of  the  results  may  lie  in  the  fact  that  the  pip  on 
the  scope  was  not  circular  but  eliptical  with  axes  of  0.2  and  0.1 
centimeters,  and  aligning  one  or  the  other  extreme  position  would 
result  in  a  d-c  error.  In  addition,  zero  set  errors  in  the  track¬ 
ing  apparatus  and  the  recording  devices  may  have  existed.  The 
foregoing  explanations  are  also  applicable  to  the  inconsistency 
with  expectations  for  £/©i  found  in  Figures  11  and  II4..  Here  one 
would  expect  the  low  frequency  values  of  the  ratio  to  be  zero,  since 
the  tracker  should  be  able  to  eliminate  the  error  in  the  d-c  case. 

Although  it  might  be  said  that  the  amplitude  ratios  for  ©0/©i, 
°oA  ,  and  e/©4  do  not  differ  significantly  over  differences  in 
subjects  and  instructions  in  view  of  the  computational  error  and 
sampling  error  confidence  limits,  the  small  amount  of  data  and  the 
few  degrees  of  freedom  available,  make  any  such  conclusions  meaning¬ 
less  at  this  time.  It  may,  however,  be  seen  that  regardless  of 
whether  there  is  any  validity  in  an  attempt  to  obtain  a  transfer 
function  from  equation  (1),  the  spectral  densities  of  output  and 
input  signals  afford  an  insight  into  the  nature  of  the  operator's 
tracking  response.  Further  research  with  different  programs,  control 
devices,  and  types  of  input  using  statistically  sound  numbers  of 
subjects  as  well  as  the  computation  of  cross  spectral  densities 
is  needed  before  this  method  can  be  properly  evaluated. 
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